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We report the first direct measurement of decays of the Af baryon involving the neutron. The 
analysis is performed using 567 pb ! of ete~ collision data collected at ys = 4.599 GeV with 
the BESIII detector at the BEPCII collider. We observe the decay At — nK$-* and measure 
the absolute branching fraction to be B(A} + nKgat) = (1.82 + 0.23(stat) + 0.11(syst))%. A 
comparison to B(A — p(K7)°) provides an important test of isospin symmetry and final state 
interactions. 
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PACS numbers: 13.30.Eg, 14.20.Lq, 13.66.Bc 


The ground-state charmed baryon A? decays eventu- 
ally into a proton or a neutron, each taking about half of 
the total branching fraction (BF) fay. However, to date 
no direct measurement of the decay modes involving a 
neutron has been performed. It has been argued that 
isospin symmetry works well in the charmed baryon sec- 
tor [2]. Comparing BFs of the final states with a neutron 
to the final states with a proton provides an important 
observable in testing isospin symmetry in At three-body 
decays [2]. The decay At — nK?z* is the most favored 
decay of the A; involving a neutron. Under the isospin 
symmetry, its amplitude is related to those of the most 
favored proton modes At + pK~at and At — pK°n°® 
as A(nK°nt) + A(pK-1+) + V2A(pK?n?) = 0. Hence, 
precise measurement of the BF for At — nK°x* provide 
stringent test on the isospin symmetry in the charmed 
baryon decays by examining this triangle relation. 


Furthermore, study of At + nK?r* is important to 
explore the decay mechanism of the A7, especially the 
factorization scheme and the involved final state interac- 
tion [21[3]. In the three-body At decay to NK, the total 
decay amplitudes can be decomposed into two isospin 
amplitudes of the NK system as isosinglet (7) and 
isospin-one (I“)). In the factorization limit, the color- 
allowed tree diagram, in which the zt is emitted and the 
NK is an isosinglet, dominates I), and I is expect- 
ed to be small compared to I as it can only proceed 
through the color-suppressed tree diagrams. Though the 
factorization scheme is spoiled in charmed meson decays, 
whether this scheme is valid in the charmed baryon At 
decays is of great interest to both theorists and experi- 
mentalists and strongly deserves the experimental inves 
tigation. The measurement of BF for At — nK?sz* 
can validate or falsify this scheme. "Together with the 
At + p(Km) the At — nK°xt can be used to de- 
termine the magnitudes of the two isospin amplitudes 
and their phase difference, which provides crucial infor- 
mation on the final state interaction. In addition, hight 
statistics data will facilitate to understand the resonant 
structures [4 [5] in the three-body A. decays and test 
the SU(3) flavor symmetry [3. Throughout the paper, 
charge conjugate modes are always implied. 


This Letter reports on the observation of the final 
states with a neutron At — nK$-^. The data ana- 
lyzed correspond to 566.93 + 0.11 pb-! [6] of e*e- an- 
nihilations accumulated with the BESIII experiment at 
V/s = 4.599 GeV H1. This energy is slightly above the 
mass threshold of a At Az pair, at which At AZ are pro- 
duced in pairs and no additional hadron is kinematical- 
ly allowed. The analysis technique in this work, which 
was first applied in the Mark III experiment [$], is spe- 
cific for charm hadron pairs produced near threshold. 
First, we select a data sample of A7 baryons by recon- 
structing exclusive hadronic decays, called the single tag 
(ST) sample. Then, we search for At — nK$m* in 


the system recoiling against the ST AZ baryons, called 
the double tag (DT) sample. In the final state nK9*, 
the neutron is not detected, and its kinematics is de- 
duced by four-momenta conservation. The absolute BF 
of At — nK&7* is then determined from the probability 
of detecting the process At — nK grt in the ST sample. 
This method provides a clean and straightforward BF 
measurement independent of the total number of At Az 
events produced. 

The BESIII detector is a cylindrical detector with a 
solid-angle coverage of 93% of 4r that operates at the 
BEPCII collider. It consists of a Helium-gas based main 
drift chamber (MDC), a plastic scintillator time-of-flight 
(TOF) system, a CsI (TI) electromagnetic calorimeter 
(EMC), a superconducting solenoid providing a 1.0T 
magnetic field and a muon counter. The charged particle 
momentum resolution is 0.5% at a transverse momen- 
tum of 1 GeV/c. The photon energy resolution in EMC 
is 2.5% in the barrel and 5.0% in the end-caps at energies 
of 1 GeV. More details about the design and performance 
of the detector are given in Ref. [9]. 

A GEANTA-based Monte Carlo (MC) simulation 
package, which includes a description of the detector ge- 
ometry and the detector response, is used to determine 
the detection efficiency and to estimate potential back- 
grounds. Signal MC samples of a A? baryon decaying 
only to nK9m* together with a A; decaying only to the 
studied tag modes are generated by the MC event gen- 
erator KKMC using EVTGEN L3. including the 
effects of initial-state radiation (ISR) Final-state 
radiation (FSR) off the charged tracks is simulated with 
the PHOTOS package [14]. The At > nK9z* decay 
is simulated using a phase space model since the two- 
body invariant mass spectra found in data for Mnr+, 
Mha K? and M Kont show no obvious structure. To study 
backgrounds, inclusive MC samples consisting of gener- 


ic A+ A> events, DoD + X production, ISR return 
to the charmonium(-like) v» states at lower masses, and 
QED processes are generated. All decay modes of the 
Ac, Y and Dis) as specified in the Particle Data Group 
(PDG) [1] are simulated by the EVTGEN MC generator, 
while the unknown decays of the w states are generated 
with LUNDCHARM [15]. 

The ST A; baryons are reconstructed using eleven 
hadronic decay modes as listed in the first column of 
Table [I] where the intermediate particles K}, A, X9, 
X and 7° are reconstructed through their decays of 
K? > ntr, À > prt, E y^ with A > prt, 
E` — pn? and 7° — yy, respectively. 

Charged tracks are required to have polar angles with- 
in | cos 0| < 0.93, where @ is the polar angle of the charged 
track with respect to the beam direction. Their distances 
of closest approach to the interaction point (IP) are re- 
quired to be less than 10 cm along the beam direction and 
less than 1cm in the perpendicular plane. Tracks origi- 
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nating from Kg and A decays are not subjected to these 
distance requirements. To discriminate pions from kaons, 
the specific ionization energy loss (dE/dx) in the MDC 
and TOF information are used to obtain particle identi- 
fication (PID) probabilities for the pion (£+) and kaon 
(£x) hypotheses. Pion and kaon candidates are selected 
using Lr > Lg and LK > Lr, respectively. For proton 
identification, information from dE/dx, TOF, and EMC 
are combined to calculate the PID probability £’, and 
a charged track satisfying £', > L'r and £', > £'x is 
identified as a proton candidate. 

Photon candidates are reconstructed from isolated 
clusters in the EMC in the regions | cos0| € 0.80 (barrel) 
and 0.86 < | cos @| < 0.92 (end cap). The deposited ener- 
gy of a neutral cluster is required to be larger than 25 (50) 
MeV in barrel(end cap) region, and the angle between the 
photon candidate and the nearest charged track must be 
larger than 10°. To suppress electronic noise and energy 
deposits unrelated to the events, the difference between 
the EMC time and the event start time is required to 
be within (0, 700) ns. To reconstruct 7° candidates, the 
invariant mass of the accepted photon pair is required to 
be within (0.110, 0.155) GeV/c”. A kinematic fit is per- 
formed to constrain the yy invariant mass to the nominal 
5? mass [I], and the x? of the kinematic fit is required to 
be less than 20. The fitted momenta of the 7° are used 
in the further analysis. 

To reconstruct K9 and A candidates, a vertex- 
constrained fit is applied to 1 * 7^ and px ^ combinations, 
and the fitted track parameters are used in the further 
analysis. The signed decay length L of the secondary 
vertex to the IP is also required to be larger than ze- 
ro. The same PID requirements as mentioned before are 
applied to the proton candidate, but not to the 7 can- 
didate. The invariant masses M,+,-, Mpr+, M,A and 
Mpo are required to be within (0.485, 0.510) GeV/c?, 
(1.110, 1.121) GeV/c?, (1.179, 1.205) GeV/c? and 
(1.173, 1.200) GeV/c? to select candidates for K}, A, 
3? and =~ candidates, respectively. 

For the ST mode pK9n°, the backgrounds involving A 
and X~ are rejected by rejecting any event with Mpr+ € 
(1.105, 1.125) GeV/c? and Mpro € (1.173, 1.200) GeV/c?. 
For the ST modes of Artnr" and $` rtr", the back- 
grounds involving KÌ and A as intermediate states are 
suppressed by requiring M,+,- € (0.480,0.520) GeV/c? 
and Mg, ¢ (1.105,1.125) GeV/c?. 


The ST A- 


c 
the variable of beam constrained mass, Mpc 


B sam 
and Pa- is the momentum of the AZ candidate. To 
improve the signal purity, the energy difference AE — 
Eveam — Ex- for each candidate is required to be with- 
in approximately +30, around the AF peak, where 
CAE is the AE resolution and Ex- is the reconstruct- 


signal candidates are identified using 
2 = 
-C = 


— Ep A-- C, where Epeam is the beam energy 


ed A; energy. The explicit AE requirements for the 
different modes are listed in Table [I] The yield of each 
tag mode is obtained from fits to the Mpc distributions 


TABLE I. ST modes, AE requirements and ST yields N;- 
in data. The errors are statistical only. 


Mode AE (GeV) Nx- 

pKZ —0.025, 0.028] 1066 + 33 
pK^«-  [-0.019, 0.023] 5692 + 88 
pKor —0.035, 0.049} 593+ 41 
pKtn «? [-0.044, 0.052] 1547+ 61 
pKonta [-0.029, 0.032] 516+ 34 
Any —0.033, 0.035} 593 + 25 
Ana —0.037, 0.052] 1864 4- 56 
Az «'«x- [—0.028, 0.030] 674-36 
von —0.029, 0.032] 532+ 30 
Er’ —0.038, 0.062] 329 + 28 
X-«*«- [—0.049, 0.054] 1009 + 57 
All tags 14415 £ 159 


in the signal region (2.280,2.296) GeV/c?, which is the 
same as in Ref. [19]. The yields of reconstructed singly 
tagged A; baryons are listed in Table [I] Finally, we ob- 
tain the total ST yield summed over all 11 modes to be 
Ne = 14415 + 159, where the error is statistical only. 


Candidates for the decay A? — nKQn* are select- 
ed from the remaining tracks recoiling against the ST 


A7 candidates. A pion with charge opposite to the ST 


Az is selected, and a K$ candidate is selected with the 
same selection criteria as described above but without 
the M,,+,- mass requirement. If more than one K9 can- 
didate is formed, the one with the largest decay length 
significance L/or, is retained, where or, is the vertex res- 
olution of L. 

Since the neutron is not detected, we use a kinematic 
variable 


Miss = Ee id |P ise] fe 
to obtain information on the missing neutron, where 
Emiss and puis; are the missing energy and momentum 
carried by the neutron, respectively, which are calculated 
by Emiss = Ebeam Exe Ec and puiss = Dat — Po = 
P+, where Ph+ is the momentum of the A? baryon, Exo 
(Pro) and Ez (Pr+) are the energies (momenta) of the 


K$ and q*, respectively. Here, the momentum p,+ is 
given by D+ = —Btagy/ Epeam/ C. — mie, where Ptag is 
the direction of the momentum of the ST A7 and ma - 


is the nominal AZ mass [i]. If the K and «+ from the 
decay A? — nKr* are correctly identified, the MZ... 
is expected to peak around the nominal neutron mass 
squared. 

The scatter plot of M,+,- versus M2... for the At > 
nKgxt candidates in data is shown in Fig. [I] where 
a cluster of events in the signal region is clearly visi- 
ble. According to MC simulations, the dominant back- 
grounds are from the decays At > X «*-«* and At > 


“tartar with M — nz+, which have the same fi- 


201711.02314v1 


chinaXiv 


0.54 |— 
«^ 0.52|- 
2 
> 
© 
e 0.5 . 
* 
R 
= 048- 
0.46 — . 
l l zi A l 
0.7 0. e 0.9 1 1.1 
PIP 
MŽ iss (GeV"/c^) 


FIG. 1. Scatter plot of M,+,- versus MŽiss for At — 
nKenxt observed from data. 


nal state as signal. Thes background events form a 
peaking background in Mĉiss, but are distributed flat in 
M-,-. Backgrounds from non-A% decays are estimated 
by examining the ST candidates in the Mpc sideband 
(2.252, 2.272) GeV/c? in data, whose area is 1.6 times 
larger than the background area in the signal region. 

To obtain the yield of At > nK$m* events, we per- 
form a two-dimensional unbinned maximum likelihood fit 
to the M2... and M,+,- distributions in both Mpc sig- 
nal and sideband regions simultaneously. As verified with 
MC simulations, we model the M,+,- and Mĉiss distri- 
butions with a product of two one-dimensional probabil- 
ity density functions, one for each dimension. The signal 
functions for M2,.. and M,+,- are both described by 
double Gaussian functions. The peaking background in 
the M2... distribution is described by a double Gaussian 
function with parameters fixed according to MC simula- 
tions, and the flat distribution in the M,+,- spectrum 
is described by a constant function. The non-A? de- 
cay background is modelled by a second-order polyno- 
mial function in the M2, distribution and a Gaussian 
function plus a second-order polynomial function in the 
M,+,- distribution, in which the parameters and the 
normalized background yields are constrained by the 
events in Mgc sideband in the simultaneous fit. The 
fit procedure is validated by analyzing a large ensemble 
of MC-simulated samples, in which the pull distribution 
of the fitted yields is in good agreement with the normal 
distribution. Projections of the final fit to data are shown 
in Fig.[2] From the fit, we obtain NP gos = 83.2 + 10.6, 
where the error is statistical only. 

The absolute branching fraction for AT 


c 
determined by 


— nK$m* is 


No Kon H 
NS X EnKont x BURG > 


B(At > nKent) = 
( c TV ST ) ata)’ 
(1) 
where En Kint İS the detection efficiency for the At — 


nK oa decay, which does not include the branching frac- 
tion for K} — ntr. For each ST mode i, the eff- 
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FIG. 2. Simultaneous fit to M2is and M,+,— of events in 
(a, b) the Az signal region and (c, d) sideband regions. Data 
are shown as the dots with error bars. The long-dashed lines 
(blue) show the AZ backgrounds while the dot-dashed curves 
(pink) show the non-A? backgrounds. The (red) solid curves 
show the total fit. The (yellow) shaded area show the MC 
simulated backgrounds from Af decay. 


ciency é K&nt is obtained by dividing the DT eliclency 
by the ST efficiency &4,. Weighting € 


i 
tag n KÜgt nK gat 


by the ST yields in data for each tag mode, we obtain 


Enkont = (45.9+0.3)%. Inserting the values of Nes 


Nx, € Enkgn+ and B(K$ — m*) [Il] in Eq. (I), we ob- 
tain B(At — nK$8s*) = (1.82 + 0.23)%, where the sta- 


tistical error, including those from NS, , and Nt is 
nKk gr As 


presented. 


With the D'T technique, the systematic uncertainties 
from the ST side cancel in the branching fraction mea- 
surement. The systematic uncertainties for measuring 
B(A} — nK&n*) mainly arise from the uncertainties 
of PID, tracking, K$ reconstruction and the fit proce- 
dure. Throughout this paragraph, all quoted system- 
atic uncertainties are relative uncertainties. The un- 
certainties in the m PID and tracking are both deter- 
mined to be 1.0% by studying a set of control sam- 
ples ofe*e"- —^mn'm ntr, ete —K*K nn and 
ete- — ppn^m«- based on data taken at energies above 
4.0 GeV. The uncertainty in the efficiency of K9 recon- 
struction is determined to be 1.596 by studying the con- 
trol samples of J/ > K** K* and J/ > óKOK*r7. 
The uncertainty due to the fit procedure is estimat- 
ed to be 5.296 by varying the fit range, the shapes of 
background and signal components, and the choice of 
sideband regions. Besides these uncertainties mentioned 
above, there are systematic uncertainties from the quot- 
ed branching fraction for Kg — mm (0.1%), the Nit 
(1.096) evaluated by using alternative signal shapes in fits 
to the Mgc spectra, the MC statistics (0.6%), the signal 
MC model (1.396) estimated by taking into account the 
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TABLE II. Summary of the relative systematic uncertainties 
for B(A + nK8n*). 


Source Uncertainty 
7 PID 1.0% 
a tracking 1.0% 
K? reconstruction 1.5% 
Fit 5.2% 
B(K3 > 1*17) 0.1% 
ve 1.0% 
MC statistics 0.6% 
MC Model 1.3% 
Total 5.9% 


statistical variations in the Mn+, Mango and Mgo;. 
spectra observed in data. These systematic uncertainties 
are summarized in Table[I]] and the total systematic er- 
ror is estimated to be 5.9% by adding up all the sources 
in quadrature. 

In summary, using 567 pb-! of e*e- collision data 
taken at ys = 4.599 GeV with the BESIII detector, 
we report the observation of the decay A? > nK$r^. 
We measure the absolute branching fraction for At — 
nKÌ%n*, B(A} + nKłr+) = (1.82+0.23+0.11)%, where 
the first uncertainty is statistical and the second is sys- 
tematic. This is the first direct measurement of a At 
decay involving the neutron in the final state since the 
discovery of the At more than 30 years ago. Quoting 
B(A] — pK «*) and B(A} — pK97°) measured by 
BESIII ia. it can be found that the amplitudes of the 
above three decay processes satisfy the triangle relation 
and validate the isospin symmetry |2]. Besides, we obtain 
B(At > nK°nt)/B(At — pK-m*) = 0.62 + 0.09 and 


B(At > nK9?z*)/B(At — pK°r®)) = 0.97 + 0.16 [18], 
in which the common uncertainties have been cancelled 
in the calculation. According to Ref. id, based on these 
ratios, the strong phase difference of J and I(9 is 
calculated to be cosó = —0.24 + 0.08, which is use- 
ful to understand the final state interactions in A? de- 
cays. Furthermore, the relative size of the two ampli- 
tudes |I(?|/|I(?| is evaluated to be 1.14 + 0.11, which 
indicates that the amplitude J“) is not small as expect- 
ed in the factorization scheme. This is consistent with 
the behaviors in the charmed meson decays [19]. These 
results will be essential inputs for the study of other Ae 
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